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Abstract. We show that the 'quantum' black hole wave operator in the k- 
Minkowski or bicrossproduct model quantum spacetime introduced in [I] has 
a natural scaling limit X p — ¥ at the event horizon. Here \ p is the Planck time 
and the geometry at the event horizon in Planck length is maintained at the 
same time as the limit is taken, resulting in a classical theory with quantum 
gravity remnants. Among the features is a frequency-dependent 'skin' of some 
— Planck lengths just inside the event horizon for lo > and just outside for 
u> < 0, where v is the frequency associated to the Schwarzschild radius. We use 
bessel and hypcrgcometric functions to analyse propagation through the event 
horizon and skin in both directions. The analysis confirms a finite redshift at 
the horizon for positive frequency modes in the exterior. 



1. Introduction 

First applications of quantum groups and noncommutative geometry to model 
quantum gravity effects have been around for some 2 decades now[2]. The work 
3J proposed how noncommutative spacetime might tame infinities in physics while 
in f 994 we introduced the particular bicrossproduct or K-Minkowski spacetime [4] 
which has since attracted a fair amount of attention, mainly because it 'almost 
commutative' and hence easier to work with. In this model the spatial generators 
Xi mutually commute but the time variable t does not. Rather, 

(1) [xi,t] = iXpX, 

where X p is a real parameter. If the noncommutative geometry arises as a model 
of quantum gravity corrections to classical geometry then presumably A p ~ 10 _44 s 
as Planck time. The model has a bicrossproduct form of quantum Poincare group 
and Section 2 provided a little more background to its construction. 

In spite of its success, which includes the proposal of its testing by time of flight 
data being collected for gamma-ray bursts by FERMI-GLAST currently in orbit, 
the main criticism of the bicrossproduct model remains that it is rather special to 
flat spacetime, relying as it does on (quantum group) Fourier transform and plane 
waves j5j. This does not sit too well as a quantum gravity effect as it effectively 
applies in some weak field limit. In [1] we have proposed an answer to this criticism; 
a general construction of a noncommutative spacetime C(M) >i T IR in which space is 
now any classical Riemannian 3- manifold (M, g) equipped with a conformal Killing 
vector field r and a static metric of classical form 

(2) /3 _1 dt(§)cr£ + g 
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where (3 £ C{M) is a function on M. In our new theory the noncommutative 
differential geometry replaces the quantum Poincare group as the guiding force, 
though the latter would be expected in the fullness of time to reemerge as the fiber 
of an associated affine quantum group frame bundle using the formalism of [5J. 
In fact the noncommutative geometry does not contain a full development of the 
noncommuative metric and noncommutative Riemannian geometry although again 
this should emerge in the fullness of time. Rather, it takes a shortcut straight to 
the physically important noncommutative spacetime wave operator. 

This theory constructs the noncommutative wave-operator on C(M ) X T K viewed 
as quantizing any static spacetime. This includes the Schwarzschild black hole. 
However, computations are quite formidable and hence in practice we work with a 
simplified version in which we use the flat bicrossproduct model spacetime ([!]) but 
'bolt on' the nontrivial Riemannian 3-geometry by a process of minimal coupling. 
This allows us to treat the black hole using the flat spacetime bicrossproduct model 
as the ambient noncommutative space. The derivation in jT] is outlined briefly in 
Section 3 to the point of the noncommutative wave equation ([6| . 

The present paper studies the black hole noncommutative wave operator in a 
certain scaling limit where distance p in Planck lengths from the event horizon is 
maintained while at the same time sending A p — > 0. This leads to a classical and 
much more amenable wave equation 

where the parameter v is the Schwarzschild frequency. This still retains quantum 
gravity effects and has the same qualitative behaviour around the horizon as the 
original theory in pQ. Our findings are summarized at the end of the paper, see 
Figure 3. The physical significance of the novel 'noncommutative skin' layer at the 
event horizon remains more fully to be understood. It is perfectly possible to read 
the paper starting at the scaling limit wave equation Q. 

2. Background to the bicrossproduct model quantum spacetime 

We first explain some of the background to the model. Indeed, quantum groups 
arose out of physics in two different classes. The first and most well known are 
the Drinfeld-Jimbo quantum groups U q (g) associated to complex semisimple Lie 
algebras g and coming out of integrable systems. Professor Jimbo at this conference 
rightly obtained the Weyl prize in part for that. The other less well-known class 
are the 'bicrossproduct' quantum groups introduced by the author as associated to 
Lie group factorisations and coming out of ideas for quantum gravity and Planck 
scale physics [2]- At the time Planck scale physics was considered unreachable by 
experiment but in recent years this view has changed and these bicrossproduct 
quantum groups have gained in interest particularly as deformations of Poincare 
and other inhomogenenous non-simple groups. 

The general construction here has semidirect product form C(M) ►<Z7(g) where 
M is a Lie group on which a Lie algebra g acts and where M acts back on g and on 
the set of its associated group G so as to form a 'matched pair'. The latter 'back- 
reaction' induces a semdirect coproduct or, equivalently, a dual quantum group 
U(m)x C(G), where m is the Lie algebra of M, making everything symmetric be- 
tween the two factors. There is also a canonical action of the first quantum group 
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on U(m) and its of dual on U(g), which are therefore naturally arising noncomma- 
tive geometries with these quantum groups as symmetries. Such data arise from 
the local factorisation of any group X into GM where neither subgroup need be 
normal, a very common occurrence in mathematics and physics. 

As well as the general construction, including specifics for the Iwasawa factori- 
sation, my 1988 PhD thesis explicitly covered [21 HE] the case 50(1, 3) « SO(3)H 3 
where H3 = M 2 x K is a curved (nonAbelian) version of M 3 . Its Lie algebra /13 
has relations ([!]) for i = 1,2, t = X3 providing a noncommutative space on which 
the bicrossproduct quantum group C(%) ^<U(so3) acts. The 4D case studied a 
bit later in [JJ is similarly based on SO(2, 3) s» £0(1, 3)7^1^ where the Lie algebra 
hi t 3 again has the relations (fTl) but with i = 1,2,3. The associated bicrossproduct 
quantum Poincare group C{Hi 3) >-<U(soi 3 ) was shown to be isomorphic as an 
abstract Hopf algebra to a 'K-Poincare quantum group' obtained by contraction of 
U q (so(2, 3)) in [2j but had a different and inequivalent interpretation of the genera- 
tors and therefore should not be confused with that. Moreover, the bicrossproduct 
model not only provided the correct interpretation of the Hopf algebra but provided 
the quantum spacetime ([lj on which it acts as the enveloping algebra U(hi^). This 
is the bicrossproduct spacetime or '/^Minkowski' model introduced by the author 
and H. Ruegg[l]. 

Since then, the bicrossproduct model quantum spacetime has attracted a great 
deal of attention, not least due to its testable variable speed of light prediction[5]. 
That there might be such a prediction had been hoped for on the basis that the de- 
formed relations of the bicrossproduct Poincare algebra involve a deformed form of 
the Casimir, but this of itself has no meaning as one can just rename the generators; 
merely calling a generator Pq or Pi does not make it energy or momentum. How- 
ever, such speculations were turned into an actual prediction in [5] by analysing the 
associated wave operator on noncommutative spacetime and its action on noncom- 
mutative plane waves e'^^e'"'. We made the reasonable assumption that physical 
plane waves would, in this model of quantum gravity effects, be identified with 
these normal ordered mathematical plane waves. 

Another feature of the bicrossproduct family is that the classical data of a 
matched pair of actions is entirely classical, consisting in the Poincare case of a 
certain non-linear action of the Lorentz group on the nonAbelian momentum group 
T-Li } 3 and vice versa. As a result the orbits of the Lorentz group in momentum 
space are classical but 'squashed' into regions with limiting surfaces or accumula- 
tion boundaries. This phenomenon is identical to that first discovered in the iso(3) 
case [3 E]j where the spherical orbits were deformed into non-concentrically 
nested spheres squashed into a region of M 3 ; in the iso(l, 3) case the mass hyper- 
bolae are similarly squashed into a cylinder |fc| < h/(c\ p ). In recent years some 
authors have sought to rebrand the bicrossproduct model as 'doubly special' be- 
cause of this remarkable feature of squashed orbits and the feature of being able 
to refer everything to classical nonlinear equations, but without actually producing 
any genuinely new model with quantum group symmetry. Both features are in fact 
typical of bicrossproduct models with noncompact factorising groups if one wants 
other examples, and I refer to [10) for a recent review. Of perhaps more interest 
is how the bicrossproduct model emerges as a limit of a q-deformed geometry and 
how this in turn emerges from quantum gravity with cosmological constant, matters 
which are explored in the 2+1 case in [IT] . 
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3. Background to the noncommutative black hole wave operator 

Given a manifold M, a metric g with Levi-Civita connection V and a conformal 
Killing vector field r, we define in pQ a noncommutative spacetime coordinate 
algebra C(M) X T K consisting of elements of the form tp = J^V'n*" where ip n G 
C{M) by which we mean, eg, smooth functions on M, and product defined via 
commutation relations 

[f,t]=t\ p T{f), V/GC(M). 

One can think of this as a cross product of C(M) by the group K with action 
defined by r and at least when M is compact make everything precise as a cross 
product C*-algebra. Next, given any pointwise invertible function f3 £ C(M) we 
define a differential calculus as, by definition, a space of 1-forms 1 (C(M) xi T M) 
spanned by usual 1-forms from M along with two extra 1-forms dt, 9', with the 
commutation relations between functions and differentials pQ 

[f,g] = 0, [f,t] = i\ P r(f), [df,g] = tX^-\df,dg)^ 

[9',f]=0, [6',t] = mX p 6', [df,t}=iX p (dr(f)-df) 
(3) [f,dt] = i\ p df, [dt,t] =tX p (0 -dt). 

for all gj £ C{M), where a = \ div(r)- 1 if M has dimension n. One of the Jacobi 
identities will quickly lead you to see that you need r to be a conformal Killing 
vector. We do rather more in PQ, building fi 1 directly on the space of classical 
1-forms f2 1 (M) extended by two extra dimensions 8',dt. The general principle of 
defining a differential structure via a bimodule of 1-forms (or an entire differential 
graded algebra) is common to most approaches to noncommutative geometry [13]. 

We also consider 2- forms 51 2 (C(M) xE) and show in [T] that the calculus is locally 
inner, i.e. in the neighbourhood of any point one can find a function h £ C(M) 
such that 

9 = dt- he', [tp, 9] = iX p d?P, {u>, 9} = -iX p duj 

for all ip in the noncommutative coordinate algebra and all noncommutative 1- 
forms cj. Finally, we show that for all normal ordered ip £ C(M) X T K the exterior 
derivative has the form 

dV> = d<A + d^dt + i^9' {A LB - ^"^^(d/W + «A p AoV^ 

where d is the classical spatial exterior derivative, do is the finite difference operator 

m - w - ix p ) 

doip{t) = . — 

i\ p 

and Alb is the classical Laplace-Betrami operator on M. We use the inverse metric 
to view the classical 1-form d/3 as a vector field, which we regard as a differential 
operator. Finally 2Aq is a deformation of the Laplacian or double-derivative in the 
t direction. We define the spacetime wave operator by 

dip = d^ + d ^dt + i^ni>9' 

in keeping with a philosophy that the extra dimension 9' plays the role of the 
conjugate to the wave operator regarded as a 'noncommutative vector field'. 
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One can apply this to the spherically symmetric 3-geometry (M, g) underlying 
the Schwarzschild black hole. Here g = h(r) 2 dr(&dr + r 2 d 2 Vt where 

a 1 2GM 

r 

for a black hole of mass M centred at the origin. Then the classical black hole has 
the form of our static metric with 

w >—w=Tr 

Moreover, the 3-geometry has a conformal killing vector field r = j^JL- Hence by 
the above we have a quantum spacetime and a wave operator on it deforming that 
of a classical black hole. 

As this wave operator and noncommutative geometry is a little hard to compute, 
we also propose £Q and in practice work with a more down to earth version of the 
theory where we let M = M 3 with its conformal Killing vector field r = r-i-. 
Then the spacetime is the usual flat spacetime bicrossproduct model with relations 
([l}. But for the calculus we take (3 as in Q which is no longer the 5D fiat space 
calculus [12] for this model. Its wave operator is then not the flat space one used 
in [5] but nor is it quite the black hole as we took a flat 3-geometry. To fix this 
we use the correct Alb ~ |/5 _1 <? _1 (d/?) in place of the one computed from the flat 
metric. We call this the 'minimally coupled black hole' as we basically adapted a 
flat space computation by replacing derivatives by covariant ones. Explicitly, 

(5) D BH m = 2A ^(t) + f( 2 - + (1 - 2) ^ + eid) i>(t + i\ p ) 

\ r or r or z J 

is our 'minimally coupled' noncommutative black hole wave operator. Here e$ = 
a x± d ^ s 2-sphere covariant derivative such that e,e ? - = — Ki+ll on the 

OXi r or L T 

spherical harmonics Y^. The nontrivial part is to compute A . Let Lip(t) — 
ip(t + i\ p ). Then[T] 

^0 = ^ (-X + ^ + (1 - 1)(X I ln(l - r ^))) , X = zA p5o . 

If we denote the expression in large parentheses as T>(X) then the noncommutative 
wave equation becomes [1] 

(2 2 7 d id 2 \ 

(6) »M) + (- - + a - + ^) V- - o 

on normal-ordered tp with ^-dependence to the right. We shall study this wave 
equation under the physical hypothesis that observed fields are to be identified in 
our model with normal ordered fields as a measure of quantum gravity corrections. 

Note that as X — > we have T>[X) = 2 (i-^) + 0(X 3 ) so we obtain — c 2^_i^ d 2 
in this limit. This has the correct classical limit for the black hole wave equation. 
Also in the limit r — !• oo we have T>(X) = — 4p + O(^) so we obtain the standard 
bicrossproduct flat space wave operator used in [S] , which replaces time derivatives 
by the finite difference ones do. 

However, the function T>(X) also has a nice limit T>(X) = — X + 4p + 0(r — 7) as 
r — > 7, i.e. close to the event horizon. So the noncommutative geometry as visible in 
the wave operator smooths out the coordinate singularity normally occurring there; 
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uj > u < 

Figure 1. Picture of the noncommutative black hole. The coordi- 
nate singularity at the event horizon is resolved into two boundaries 
with a frequency dependent 'interregnum' or skin between them of 
depth 7|e~" Ap — 1| where 7 is the Schwarzschild radius. The picture 
is different for positive and negative frequency uj modes. 



in some sense it makes it one dimension better as a log in place of a l/(r — 7). It 
similarly improves matters at the other boundary r = je~ uXp which is the other side 
if an interregnum or 'skin' region in which the logarithm has a negative argument. 
We will take a close look at both boundaries by means of a scaling limit. The 
geometrical picture that we refer to is shown in Figure 1 for positive and negative 

UJ. 



4. Solutions of the wave equation near the classical event horizon 
To study this limit r — > 7 further, we set 

r-j = l p p 

where l p — c\ p is the Planck length and p is the distance from the horizon in Planck 
units. Then 

- i[ —)d . 



1-? C P 
Hence rl p times the wave equation ^ becomes 



l p cp 7 + l p p dp dp 2 

Although lp is what it is (about 10 -33 cra) we can simplify this equation from a 
mathematical point of view by setting l p — ¥ 0. This should be viewed as no more 
than a certain scaling limit approximation in which some quantum gravity effects 
are ignored but a certain remnant remains. In this limit the wave equation in these 
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( 8 ) ( -2 ( -- + pln(l + — )) + jr.+Pj^ )^ = 



coordinates becomes 

if we assume bounded angular momentum so that — > 0. Here 

c 

^ = — 

7 

is, roughly speaking, the frequency of a wave of wavelength the size of the black 
hole, as 7 is the Schwarzschild radius. 

We look at tp with time-dependence e lult (separating variables). Then the equa- 
tion Q becomes for each mode of frequency ui, 

d_ 

v r v pv' J ' dp ^ dp 2 

where ip is now a function just of p. Remembering the original physical interpre- 
tation of p, we see that for a given sign of w there is a strip of width \u)\/v Planck 
lengths where the argument of the logarithm is negative. Specifically, this happens 
if 

M M 

uj>0: -^<p<0, oj<0: < p < — . 

V V 

Thus we have qualitatively the same kind of 'quantum skin' or 'interregnum' region 
in our scaling limit as for the original theory in Figure 1, but now with width \ui\jv. 
This lies just inside the event horizon for positive frequencies and just outside it for 
negative frequencies. As long as we stay out of these regions, we have real solutions. 
We use MATHEMATICA to compute numerical solutions, but our goal is to gain 
an analytic understanding of them in different regimes. We focus on uj > but 
return to the negative frequency case at the end. 

(i) p » " Far outside the event horizon If we assume that 

(9) M « P 

then we can expand the logarithm. At lowest nontrivial order we obtain the con- 
ventional wave equation in unusual coordinates 

_L_^!_ + — + ,1,-0 
pv 2 dt 2 dp ^ dp 2 ) 

Thus, if we write p = then the wave equation has solutions 

ip±{t,p) =e lul{t± ^ =e lut p ±l *. 

If we do a little better and expand to the next cubic order, the wave equation takes 
the approximate form 

( 1 d 2 n 1 d 3 d d 2 \ , n 

(10) \-pV 2 W 2 - 2 Jp^W + d P + p di 2 ) i) = 

This has solutions 



where I n (z) is a bessel I function but extended to complex n and the T function 
is a normalisation. In another context one might consider taking ilu = uj' real in 
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FIGURE 2. Scaling limit solutions against distance p Planck lengths 
from event horizon, (a) Numerical solution for ^ = 10 with differing 
boundary conditions from p — 10 5 to p — 10 -3 on log scale (b) Com- 
parison with real and imaginary parts bessel solution for p >> ^ = 10 
(right) and 10 < p << ^ (left), (c) Real and imaginary parts numer- 
ical solution (left) and comparison with hypergeometric solution (right) 
for ^ = 5.5 in 'skin' — ^ < p < driven by real boundary condition at 
the horizon, on linear scale, (d) Ditto numerical solution but for 77 = 1 
driven by boundary condition at the interior boundary, (e) Numeri- 
cal solution (left) for ^ — 10 with differing boundary conditions from 
p = — 10 1 to p — — 10 3 , and comparison (right) with \p\ >> — bessel 
solution, on log scale. 
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which case note the poles in the normalisation factor of ip± at =F positive integer 
=^-. Figure 2(b) on the right confirms, however, that the case of real u> is wave-like 
and indeed asymptotes to the plane waves just found when p is large. Note that 
the plot is against log p with ~ = 10. The approximation (|9j) means that we can 
only approach the horizon to uj/v Planck lengths before quantum gravity effects 
render the weak field (i.e. small 7 or large v) expansion invalid. Thus, the right 
side of Figure 2(b) is only valid for p >> 10 1 and this is confirmed on comparison 
with the numerical solution in Figure 2(a). 

(ii)0 < p << " Just outside the event horizon The numerical solution in Fig- 
ure 2(a) continues down to p = 10~ 3 or so on the log plot. At some point we enter 
the regime 

(11) \p\ « - 

v 

where the logarithm is suppressed and the wave equation becomes 

(12) 2 ^ = {wp +p wr 

As a result the wave stops being so exponentially stretched in the direction of 
increasing p (it begins to be stretched in the logarithmic plot) and starts to behave 
as a linear combination of 




where Jq,Y$ are bessel J and Y functions. These have respectively a finite value 
and zero slope at the event horizon p = or an infinite value and finite slope at 
p = in the logarithmic plot. Thus we conclude that the equation ^ has a full 
set of solutions outside the event horizon becoming exponentially stretched plane 
waves for large p and tending to either a constant or a constant p-§^ip (the slope in 
the logarithmic plot) at the event horizon p — 0. A generic boundary condition set 
at large p will stimulate the latter Jg mode and result in a log divergence at the 
horizon, so this is presumably the fate of an infalling wave. However, a wave created 
at the event horizon (which would be relevant to Beckenstein-Hawking radiation) 
could take the form of the regular Jo mode. 

(iii)— " < p < In the skin layer We now look inside the event horizon. Initially, 
for small negative p we are in the same regime as in (ii), except that there is now 
a term —2mp term in the wave operator if we take the standard choice of branch 
of the logarithm. This system can be solved using hypergeometric functions as 

* . e,* lFl (I + ,, - 2p ^) , e ^. Fu (I + ,, - 

where Fjj is the confluent hypergeometric function. We focus for the moment on \F\ 
which has a finite value and zero slope at the event horizon and can therefore match 
the Jo (bessel J) on the other side. Its continuation further into the 'skin' region 
acquires an imaginary part and both parts grow rapidly as we go deeper into the skin 
(more negative p) . The values at p = — ^ are (according to MATHEMATIC A) given 
by a Laguerre L function and appear to grow very rapidly with ^. In Figure 2(c) 
we show this analytic approximate solution on the right and for comparison the full 
numerical solution, which is similar. The boundary condition is a real value 1 and 
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zero slope at or close to p — and we have reduced to ^ — 5.5 in order to have a 
visible comparison at a point where significant 'amplification' starts to appear. 

For solutions driven by real conditions at (or rather, just inside) the interior 
boundary p — — ^ the numerical solution is shown in Figure 2(d) for a modest 
jf = 1. In general a significant visible amplification in the numerical solution (say 
a factor of some 10 6 ) sets in around ^ = 4 and a comparable imaginary part 
is similarly acquired, i.e. much as we had going the other way. However, close 
inspection of the numerical solution appears to indicate that even for the example 
shown, there is in fact a log divergence at p = 0. (This is more visible in the middle 
scenario in our later Figure 3 where we have used a log scale so that a log divergence 
is a straight line and where we have suppressed a very large imaginary component 
in the plot allowing us to see a bit more detail than in Figure 2(d).) This brings 
us back to the hypergeometric mode Fjj. This indeed has a log divergence in its 
real and imaginary part and a similar form in the region — ^ << p < 0. Even 
though the approximation used breaks down as we approach p = — — , as confirmed 
by the detail in Figure 3 mentioned, a generic boundary condition there can be 
expected to propagate into the region where the approximation ( |11[ ) is valid and 
thereafter drive this mode to produce a log divergence. The divergence can be 
matched in real slope at p = on the log plot to that of a Yq (bessel Y) mode 
outside, although the role of the similarly divergent imaginary part is not entirely 
clear. It is presumably related to the logarithmic asymptotic form. Overall, we see 
that there is a significant asymmetry in the modes crossing in the 'skin layer' from 
the black hole interior, which appear to diverge as Fjj just inside the horizon, and 
Jo modes crossing from outside the event horizon which are regular at the horizon 
and match to \F\ modes inside. 

(iv),o << — ^ Black hole interior As we continue further with negative p we 



return again to the regime \p\ » ^ and we are back to (10), this time with p < 0. 
If we change variables to — p we have the same equation but with —ui. Hence 
the cubic correction enters with the opposite sign. We have the same solutions 
Bessel functions /„ with imaginary n as before but now an i in the argument. As 
Figure 2(e) on the right confirms, we still have wave-like behaviour and these tend 
for large negative p to plane waves in log plot. The difference is that whereas the 
bessel function in regime (i) showed a slight stretching of the waves in the log plot 
compared to large p, now we see a slight compression in the log plot compared 
to large negative p. We again see a good fit with the numerical solutions, shown 
for real boundary conditions of finite amplitude/zero derivative, or vice- versa. We 
have returned to — = 10 and have matched normalisation at large negative p. 

(v) Negative frequencies We have focussed on u > 0. If we suppose instead that 
u) < then note in ^ that we can change back to u) > if we also change p to 
—p. Hence the story for uj < is exactly the same as above but with a left-right 
reflection in all of Figure 2. This time a mode entering from far outside the event 
horizon is slightly (but finitely) compressed compared to the overall frequency- 
independent compression expressed in the log plot, and diverges just above the 
classical event horizon. Meanwhile, a mode entering from the black hole 'interior' 
is slightly (but finitely) stretched in the log plot, crosses the classical event horizon 
without coordinate singularity and is amplified as it passes out the other side. 

All of this behaviour is remarkably similar to that of the full theory in [T], sug- 
gesting that our simplified scaling limit wave equation ^ shares many qualitative 
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Figure 3. Pastice of black hole wave solutions in scaling limit 
at the classical event horizon p = for u> > on a log plot. 
The system is driven by the location of the boundary conditions 
as indicated by the arrow. In the interior and exterior regions 
we have set ^ = 10 but have artificially reduced this in the skin 
regions to ^ = 1.6 in order to have a viable diagram with visible 
detail, and also note that solutions in the skin region acquire a 
comparable imaginary part which we omit for clarity. At top, 
we show a typical mode set far in the exterior, resulting in a log 
divergence but finite redshift just outside the horizon p = 0. In 
middle we show a mode set in the 'far' interior, resulting in a log 
divergence just inside the horizon. At bottom we show solutions 
set at the horizon. The pictures for u) < are the same with 
'exterior' and 'interior' interchanged, but note that the boundary 
condition illustrated would still be driving the skin region from the 
classical event horizon which is now just below the skin of the black 
hole. 
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features with the full wave equation. The only thing missing is that we cannot 
investigate 'standing waves' in the black hole interior. This is because our scal- 
ing limit covers only the geometry near the horizon and a fixed number of Planck 
lengths from it (this includes the skin region). In effect, on entering the black hole 
interior the wave solutions in Figure 2(e) continue on towards the origin but this is 
considered as infinitely far away for a macroscopic black hole. On the other hand, 
one could consider a different scaling limit in which 7 = nl p and n is fixed as l p — > 
to analyse Planck size black holes in the model. This will be looked at elsewhere. 

Meanwhile, the propagation in the two directions across the event horizon is 
summarised in Figure 3. Because the orbital angular momentum was washed out 
in our scaling limit and because we are far from the centre of the black hole, the 
system is more like propagation in a line. More precisely, two lines laid end to end. 
Thus, one could change variables to x = ^\ogp for p > and y = — 7log(— p) for 
p < and arrive at a more conventional -j^ or -j^ for the spatial part. We have 
refrained from this only in order to keep the original context. In the log plot in 
Figure 3 the interior and exterior propagations each span a factor of 10 6 in p and are 
set at ^ = 10. However, this would result in extremely high 'amplification' factors 
in the skin region and hence for purposes of illustration only, we have set " — 
1.6 for solution in the skin region. The boundary conditions are generally mixed 
ones setting amplitude and gradient, with the arrows indicating the location of the 
boundary condition and solving the equations from there, rather than indicating 
'direction of travel', which we have not particularly analysed. Although it does not 
appear to be too singular the equation at p — — ~ is numerically ill-conditioned 
and as a result the boundary conditions are set just inside the skin region rather 
than right on the boundary. 

We should also say that the physical interpretation of all of these features remains 
problematic. We have not discussed here issues such as unitarity (this is likely to 
be a problem) but it should be remembered that the entire noncommutative model 
is intended as an effective description of quantum gravity effects rather than a self- 
contained closed system. However, at least at first sight it would appear that waves 
from infinity are likely to log diverge at the classical event horizon p = 0. We have 
seen that one can also create modes at the event horizon which are regular there 
and extend to waves outside. Inside the skin they grow rapidly and acquire an 
imaginary component, and the (now complex) solution appears to be matchable 
to waves continuing into the interior. This is shown in the lower part of Figure 3. 
Meanwhile, in the middle of Figure 3 we show waves in the interior which appear to 
be matchable to a solution inside which becomes complex and has a log divergence 
at p = 0, this time on the inside of the classical event horizon. All of this is for 
positive frequency with a reversed picture for negative frequency. This represents a 
significant time-asymmetry, which is to be expected. Whether the log divergences 
at p = arc important is not clear as they could be artefacts of the scaling limit; 
in the actual quantum gravity model presumably other physical effects would take 
over. 

References 

[1] Majid S 2010, Almost commutative Riemannian geometry, I: wave operators, arXiv:1009.2201 
(math.QA), 50pp. 

[2] Majid S 1988, Hopf algebras for physics at the Planck scale, J. Class. Quant. Gravity 5, 
1587-1607. 



SCALING LIMIT OF THE NONCOMMUTATIVE BLACK HOLE 



13 



[3] Majid S 1990, On q-rcgularization, Int. J. Modern Physics A. 5, 4689-4696 

[4] Majid S& Rucgg H 1994, Bicrossproduct structure of the K-Poincarc group and non- 
commutative geometry, Phys. Lett. B. 334, 348-354 

[5] Amelino-Camelia G & Majid S 2000, Waves on noncommutativc spacctimc and gamma-ray 
bursts, Int. J. Mod. Phys. A15, 4301-4323. 

[6] Majid S 1999, Quantum and braided group Riemannian geometry, J. Geom. Phys. 30, 113-146 

[7] Majid S 1990, Matched pairs of Lie groups associated to solutions of the Yang-Baxter equa- 
tions, Pacific J. Math 141, 311-332. 

[8] Majid S 19991, Hopf-von Neumann algebra bicrossproducts, Kac algebra bicrossproducts, 
and classical Yang-Baxter equations, J. Functional Analysis 95, 291-319. 

[9] Lukierski J, Nowicki A, Ruegg H, & Tolstoy VN 1991, (jr-Dcformation of Poincarc algebra. 
Phys. Lett. B, 268, 331-338. 
[10] Majid S 2009, Algebraic approach to quantum gravity II: noncommutativc spacetimes, in 

Approaches to Quantum Gravity, ed. D. Oriti. C.U.P. 466-492 
[11] Majid S & Schroers B 2009, q-Deformation and semidualisation in 3D quantum gravity, J. 

Phys A 42, 425402 (40pp) 
[12] Sitarz A 1995, Noncommutativc differential calculus on the K-Minkowski space, Phys. Lett. B 
349, 42-48. 

[13] Connes A 1994, Noncommutative Geometry, Academic Press. 

Queen Mary University of London, School of Mathematics, Mile End Rd, London El 
4NS, UK 

E-mail address: s.majid@qmul.ac.uk 



